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Abstract 

It is well-known how to determine the price of perpetual American 
options if the underlying stock price is a time-homogeneous diffusion. 
In the present paper we consider the inverse problem, i.e. given prices 
of perpetual American options for different strikes we show how to con- 
struct a time-homogeneous model for the stock price which reproduces 
the given option prices. 



1 Introduction 

In the classical Black-Scholes model, there is a one-to-one correspondence 
between the price of an option and the volatility of the underlying stock. 
If the volatility a is assumed to be given, for example by estimation from 
historical data, then the arbitrage free option price can be calculated using 
the Black-Scholes formula. Conversely, if an option price is given, then the 
implied volatility can be obtained as the unique a that would produce this 
option price if inserted in the Black-Scholes formula. It has been well doc- 
umented that if the implied volatility is inferred from real market data for 
option prices with the same maturity date but with different strike prices, 
then typically a non-constant implied volatility is obtained. Since the im- 
plied volatility often resembles a smile if plotted against the strike price, this 
phenomenon is referred to as the smile effect. The smile effect is one indica- 
tion that the Black-Scholes assumption of normally distributed log-returns 
is too simplistic. 
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A wealth of different stock price models lias been proposed in order to 
overcome the shortcomings of the standard Black-Scholes model, of which 
the most popular are jump models and stochastic volatility models. Given a 
model, option prices can be determined as risk neutral expectations. How- 
ever, models are typically governed by a small number of parameters and 
only in exceptional circumstances can they be calibrated to perfectly fit the 
full range of options data. 

Instead, there is a growing literature which tries to reverse the procedure 
and uses option prices to make inferences about the underlying price process. 
At one extreme models exist which take a price surface as the initial value 
of a Markov process on a space of functions. In this way the Heath- Jarrow- 
Morton |Ej interest rate models can be made to perfectly fit an initial term 
structure. Such ideas inspired Dupire [5] to introduce the local volatility 
model which calibrates perfectly to an initial volatility surface. For a local 
volatility model, Dupire derived the PDE 

PT{T,K) + rKPK{T,K) = ^a\T, K)K^ Pkk{T, K) 

where P{T, K) is the European put option price, T is time to maturity and 
K is the strike price. Solving for the (unknown) local volatility cr(T, K) gives 
a formula for the time inhomogeneous local volatility in terms of derivatives 
of the observed European put option prices. 

The local volatility model gives the unique martingale diffusion which 
is consistent with observed call prices. (Alternative, non-diffusion models 
also exist, see for example, Madan and Yor [9].) The recent literature (eg 
Schweizer and Wissel [13]) has included attempts to extend the theory to 
allow for a stochastic local volatility surface. However, it relies on the knowl- 
edge of a double continuum of option prices, which are smooth. In contrast, 
Hobson [7j builds models which are consistent with a continuum of strikes, 
but at a single maturity, in which case there is no uniqueness. 

In the current article we present a method to recover a time-homogeneous 
local volatility function from perpetual American option prices. More pre- 
cisely, we assume that perpetual put option prices are observed for all dif- 
ferent values of the strike price, and we derive a time-homogeneous stock 
price process for which theoretical option prices coincide with the observed 
ones. 

No-arbitrage enforces some fundamental convexity and monotonicity 
conditions on the put prices, and if these fail then no model can support the 
observed prices. If the observed put prices are smooth then we can use the 
theory of differential equations to determine a diffusion process for which 
the theoretical perpetual put prices agree with the observed prices, and our 
key contribution in this case is to give an expression for the diffusion co- 
efficient of the underlying model in terms of the put prices. (It turns out 
that this expression uniquely determines the volatility co-efficient at price 
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levels below the current stock price, but the volatility function is undeter- 
mined above the current stock price level, except through a single integral 
condition.) The key idea is to construct a dual function to the perpetual 
put price, and then the diffusion co-efficient can be found easily by taking 
derivatives of this dual. 

The second contribution of this paper is to give time-homogeneous mod- 
els which are consistent with a given set of perpetual put prices, even when 
those put prices are not twice differentiable or not strictly convex (in the 
continuation region where it is not optimal to exercise immediately). Again 
the key is the dual function, coupled with a change of scale and a time- 
change. We give a construction of a time- homogeneous process consistent 
with put prices, which we assume to satisfy the no-arbitrage conditions, but 
otherwise has no regularity properties. 

One should perhaps note that in reality put prices are only given in 
the market for a discrete set of strike prices. Therefore, as a first step one 
needs to interpolate between the strikes. If a stock price is modeled as 
the solution to a stochastic differential equation with a continuous volatility 
function, then the perpetual put price exhibits certain regularity properties 
with respect to the strike price. Therefore, if one aims at recovering a con- 
tinuous volatility, then one has to use an interpolation method that produces 
option prices exhibiting this regularity. On the other hand, if a linear spline 
method is used, then a continuous volatility cannot be recovered. This is 
one of the motivations for searching for price processes which are consistent 
with a general perpetual put price function (which is convex but may be 
neither strictly convex, nor smooth). 

Whilst preparing this manuscript we came across a preprint by Alfonsi 
and Jourdain [2]. One of the aims of pj, as in this article, is to construct 
a time-homogeneous process which is consistent with observed put prices. 
However the method is different, and considerably less direct. Alfonsi and 
Jourdain [2j build on their previous work [1] to construct a parallel model 
such that the put price function in the original model (expressed as a func- 
tion of strike) becomes a call price function expressed as a function of the 
initial value of the stock. They then solve the perpetual pricing problem 
for this parallel model, and subject to solving a differential equation for the 
optimal exercise boundaries in this model, give an analytic formula for the 
volatility coefficient. In contrast, the approach in this paper is much sim- 
pler, and unlike the method of Alfonsi and Jourdain extends to the irregular 
case. 
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2 The Forward Problem 



Assume that the stock price process X is modeled under the pricing measure 
as the solution to the stochastic differential equation 

dXt = rXt dt + a{Xt)Xt dWt, Xq = xq. 

Here the interest rate r is a positive constant, the level-dependent volatility 
a : (0, oo) (0, oo) is a given continuous function, and W is a standard 
Brownian motion. We assume that the stock pays no dividends, and we let 
zero be an absorbing barrier for X. If the current stock price is xq, then the 
price of a perpetual put option with strike price K > is 

P{K) = supE^« [e-'^^iK - X,)+] , (1) 

where the supremum is taken over random times r that are stopping times 
with respect to the filtration generated by W . From the boundedness, mono- 
tonicity and convexity of the payoff we have: 

Proposition 2.1. The function P : [0, oo) — > [0,oo) satisfies 

(i) {K - xo)+ < P{K) < K for all K. 

(a) P is non- decreasing and convex on [0,oo). 

Example. If a is constant, i.e. if X is a geometric Brownian motion, then 

P(K) = / &(W^o(/3 + 1))^ iiK<k 

\ K -XQ ifK>k, ^ ^ 

where /3 = 2r/a'^ and K = xo{P + 1)//?. 

Intimately connected with the solution of the optimal stopping problem 
([l]) is the ordinary differential equation 

-cr{x)'^x'^Uxx + rxux — ru = 0. (3) 

This equation has two linearly independent positive solutions which are 
uniquely determined (up to multiplication with positive constants) if one 
requires one of them to be increasing and the other decreasing, see for ex- 
ample Borodin and Salminen [H pl8]. We denote the increasing solution by 
tp and the decreasing one by (p. In the current setting, ip and ip are given by 

ip{x) = X 

and 

f°° I ( n 2r } 
m = xj^ ^exp|-y^ ^dzjdy (4) 
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for some arbitrary constant c. For simplicity, and without loss of generality, 
we choose c so that (^(xq) = 1- It is well-known that with = inf{t > : 

Xt = z], 

^ \ xp{x)/i/;{z) li z > X. ^ ^ 

(This result is easy to check by considering e~'^'^0{Xt) and e~^^'ip{Xt) which, 
since they involve solutions to ([3]), are local martingales.) 

Given the assumed time-homogeneity of the process X, it is natural to 
consider stopping times in ([T]) that are hitting times. Define 

P{K) := sup E^o [e-"^^(i^-XHj+] 

z:z<xoAK 

sup {K - z)E^" [e-^^-] 

z:z<X()AK 

K-z 

= sup , 

z:z<xoAK ^{Z) 

where the last equality follows from ([5]). By differentiating we find 
= exp {- ^ d.] £ 1 (-P {- £ 7^ - l) 

(6) 

SO ip'{xo) < 0. Define 

K = xq- l/(p'{xo). 

Then, for K < K we have 

K — z K — z 

P{K)= sup — — — =sup— — — , (7) 

z:z<K V{z) z ^{Z) 

whereas for K > K we have P{K) = {K — xq). 

Clearly, P{K) > P{K), and of course, as we show below, there is equal- 
ity. 

Lemma 2.2. The functions P and P coincide, i.e. 

P{K) = sup (8) 

z:z<XQ ^\Z} 

Proof. As noted above, clearly P > P since the supremum over all stopping 
times is at least as large as the supremum over first hitting times. 

For the reverse implication, suppose first that K < K. In that case 
Lp{z) > {K — z)^ /P{K) by ([71). Further, e~'^^<p{Xt) is a non-negative local 
martingale and hence a supermartingale. Thus, for any stopping time r we 
have 

1 > E^o[e-''^(^(X^)] > E^«[e-''^(i^ - Xr)+/P{K)]. 
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Hence P{K) > sup^ E^o [e-''^(i^ - = P{K). 

Finally, let K > K. It follows from the first part that P{K) = K — xq, 
so Proposition 12.11 implies that P{K) = K — xq = P{K), which finishes the 
proof. □ 



Example. If a is constant, i.e. if X is a geometric Brownian motion, then 

X ) 

where (3 = 2r/o"^. Consequently, the put option price is given by 



P{K) = Xq^ sup {K - z)z^. 

z:z<xo 

Straightforward differentiation shows that the supremum is attained for 

if z* < Xq, and for z = xq if z* > xq. Consequently P{K) is given by ([2]). 

Proposition 2.3. There exists a K £ {0, oo) such that P{K) > {K — xo)~^ 
for all K £ (0, K) and P{K) = K - xq for all K>k. 

Proof. Since a is continuous and positive, the function 0[x) given by (jl]) is 
twice continuously differentiable with a strictly positive second derivative. 
Therefore, for each fixed K there exists a unique z = z{K) < xq for which 
the supremum in Lemma 12.21 is attained, i.e. 

PiK) = ^-^(^^ (9) 

Geometrically, z = z{K) is the unique value which makes the (negative) 
slope of the line through (i^, 0) and {z,Lp{z)) as large as possible. Alge- 
braically this means that z = z{K) and K satisfy z{K) = xq li K > K 
where 

k = xo- l/(^'(xo), 

and 

(K - z)^'(z) + (^(z) = (10) 
UK < k. Therefore, ifK>k, then 

P{K) = ^-^=K-xo, 

and ii K < k, then 

P{K) = ^^^>K-xo. 

□ 
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Proposition 2.4. Suppose a : (0, oo) (0, oo) is continuous. In addition 
to the properties in Propositions \2.1\ and \2.!^ the following statements about 
the function P : [0, oo) [0, oo) hold. 

(i) P is continuously differentiable on (0,00), and twice continuously dif- 
ferentiable on (0, 00) \ {K}. 

(a) P is strictly increasing on (0, 00) with a strictly positive second deriva- 
tive on (0, K). 

Proof. It follows from pU|) and the implicit function theorem that z{K) is 
continuously differentiable for K < K. Therefore, differentiating ([9]) gives 

^ (1 - z'{K))mK)) -{K- z{K))z'{K)^'{z{K)) 
1 



^{z{K)y 

where the second equality follows from (jlOp . Equation (jlip shows that 
P'{K—) = l/(^(xo) = 1, so P is also at K. Moreover, since ip{z) is 
and z(K) is away from K, it follows that P{K) is on (0, 00) \ {K}. 
In fact, for K < K we have 

-z'{K)ip'iz{K)) 



P"(K) = 

^ ' m^K)))^ 



{K-z{K))(^{z{KW^"{z{K)) 

where the second equality follows by differentiating pup . Thus P has a 
strictly positive second derivative on (0,i^), which finishes the proof. □ 

Remark. Note that P'(0+) > with equality if and only if <^(0+) = 00. 

We end this section by showing that <p can be recovered directly from the 
put option prices P{K), at least on the domain (0,xo]. To do this, define 
the function 93 : (0, xq] — > (0, 00) by 

K — z 

ip{z) = sup -f——, (12) 

K:K>z P{K) 

where P is given by ([8]). 

Lemma 2.5. (i) Suppose f : {0,zq] [l,oo] is a non-negative, decreas- 
ing convex function on {0,zo] with f{zo) = 1 and f'{zo) < 0. Define 
g : (0,00) [0,00) by 

g{k) = sup ^jT^-- (13) 

z:z<zo J\^) 

Then g{k) is a non-negative, non- decreasing convex function with (k — 
zo)^ < g{k) < k and g{k) = k — zq for k > k* = zq — 1/ f'{zo). 
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(ii) f and g are self-dual in the sense that if for z < zq we define 

k — z 



F{z) = sup 

k:k>z 9\f^j 



then F = f on (0, zq]. 



(Hi) Similarly, assume that g : (0, cxd) [0, oo) is a non-negative, non- 
decreasing convex function with (A; — zq)^ < g{h) < k and g{k) = k — XQ 
for k > k* , and define 

f{z) = sup ^--^ 

k:k>z 9[k) 

for z < zq. Then g can be recovered by 

g{k) = sup ^— ^. 



Proof. See Appendix A.l. □ 

Corollary 2.6. The function ip coincides with the decreasing fundamental 
solution on (0,a;o]. 



3 The Inverse Problem: The Regular Case 

Now consider the inverse problem. Let P{K) be observed perpetual put 
prices for all nonnegative values of the strike K. The idea is that since 
(p satisfies the Black-Scholes equation ([3]), Corollary 12.61 provides a way to 
recover the volatility a{x) for x G (0, xq] from perpetual put prices. In 
this section we provide the details in the case when the observed put prices 
are sufficiently regular. We assume that the observed put option price P : 
[0,oo) — > [0, oo) satisfies the following conditions, compare Propositions 12 . 1 1 
O and [Za above. 

Hypothesis 3.1. (i) {K - xo)+ < P{K) < K for all K. 

(ii) There exists a strike price K* such that P{K) > (K — xq)^ for all 
K <K* and P{K) = K - xq for all K > K* . 

(Hi) P is continuously differentiable on (0, oo), and twice continuously dif- 
ferentiable on (0, oo) \ {K*}. 

(iv) P is strictly increasing on (0, oo) with a strictly positive second deriva- 
tive on (0, K*). 
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Motivated by Corollary 12.61 define the function (p : (0,j;o] (0,cxd) by 

K — z 

ip{z) = sup -^TTTT. (14) 

K:K>z j 

Proposition 3.2. The function P can he recovered from ip by 

K-z 



P{K) = sup 

Proof. This is a consequence of part (iii) of Lemma |2.5[ □ 

Proposition 3.3. The function ip : (0, xq] — > (0, oo) is twice continuously 
differ entiable with a positive second derivative, and it satisfies 'p{xq) = 1 
and ip'ixo) = -1/{K* - xq). 

Proof. For each z < xq there exists a unique K = K{z) G {z,K*] for which 
the supremum in (jl4p is attained. Geometrically, K is the unique value 
which minimises the slope of the line through (z, 0) and {K, P{K)). Clearly, 
K = K{z) satisfies the relation 

{K - z)P'{K) = P{K). (15) 

Reasoning as in the proof of Proposition 12.41 one finds that K{z) is contin- 
uously differentiable on (0,a;o], with 

and p'{xq) = —1/{K* — xq). Differentiating fl6]) with respect to z gives 

K'{z)P'{K{z)) {P'{K{z))f 
^ ' P^K{z)) {K{z)- z)P'^{K{z))P"{K{z)y ^ ' 

where the second equality follows by differentiating (|15p . It follows that 
(p"{z) is continuous and positive, which finishes the proof. □ 

Next, define cj(x)^ for x < xq so that p is a solution to the corresponding 
Black-Scholes equation, i.e. 

Moreover, define o"(x)^ for x > xq to be the constant 

= 2r{K* -xo)/xo. (19) 

Now, given this volatility function ct(-), we are in the situation of Section [21 
and can thus define p to be the decreasing fundamental solution to the 
corresponding Black-Scholes equation scaled so that p{xq) = 1. Moreover, 
let P{K) be the corresponding perpetual put option price as given by dS]). 
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Theorem 3.4. Assume that Hypothesis \3.1\ holds. Then the functions P and 
P coincide. Consequently, the volatility cr{x) defined by ()18p - (fT9|) solves the 
inverse problem. 

Proof. First recall that any solution to ([3]) can be written as a linear com- 
bination of tp and (p. Consequently, it follows from the definition of a that 

ifix) = Ci^{x) + C2X (20) 

for X £ (0,xo] and some constants Ci and C2. The boundary condition 
if{xQ) = 1 yields 

Ci + C2X0 = 1. (21) 

Moreover, from ()19|) we have 

XQ [ ^expj- / dz\ dy 



Xo I ^ exp I - / —TT^ T dz } dy 



I Jxo - Xo) 

K*/(K*-xo) / 2-xo/(K*-xo) 



JXQ 



SO the requirement (p{xo) = 1 gives the explicit expression 
for the decreasing fundamental solution. Consequently, 



j{* / r<xi I ( fy 2r 



K* fK*-xo 1 



K* — Xo V K*xo 



-\ =-l/{K* -xo) = ip'{xo) 
xoj 



by Proposition 13.31 It therefore follows from (f20|) that (1 — Ci)(^'(2;o) = C2. 
Since (^'(2^0) < 0, (PT|) gives that Ci = 1 and C2 = 0, so = on (0,j;o]- 
Proposition 13.21 then yields 

P{K) = sup = sup = PiK), 

z:z<xo f{Z) z:z<xo f\Z) 

which finishes the proof. □ 

Remark. The inverse problem does not have a unique solution. Indeed, a 
can be defined arbitrarily for x > xq as long as the condition 

f°° 1 { n 2r \ 

/ — exp<^-/ — —-2dz\dy = {K*-xo)/{K*xo) 

Jxo y I Jxo ZCT\Z) ) 

is satisfied, which guarantees that is at the pointxQ. 
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We next show how to calculate the volatility that solves the inverse 
problem directly from the observed option prices P{K). To do that, note 
that for each fixed z, the supremum in (I14p is attained at some K = K{z) 
for which 

, ^ K — z , ^ 

^U) = (22) 

= (23) 

and 

^ '~ {K-z)P^{K)P"{Ky ^ > 

compare (fT6|) and p!7j) . Since ip satisfies the Black-Scholes equation, we get 

. _ M^)-z^'{z) _ 2rKP\K)P"{K) 

Consequently, to solve the inverse problem one first determines z by 

z = K-^ 
P'iKY 

and then for this z one determines cr(z) from (|25p . 



4 The inverse problem: the irregular case 

Again, suppose we are given perpetual put prices P{K) and an interest 
rate r > 0. Our goal is to construct a time-homogeneous process which 
is consistent with the given prices. Unlike in the regular case discussed 
in Section [3l we now impose no regularity assumptions on the function P, 
beyond the necessary conditions stated in Propositions 12.11 and 12.31 For a 
discussion of the necessity of the condition in Proposition l2.3l see Section [^TTl 

Hypothesis 4.1. (i) For all K we have {K - xo)+ < P{K) < K. 

(ii) P is non- decreasing and convex. 

(Hi) There exists K* £ (xq, oo) such that P{K) = K — xq for K > K* . 

Given P{K) we define K £ [0,xo] by ^ = sup{i^ : P{K) = 0}. Then, 
necessarily, P{K) = for K < K^. 

Theorem 4.2. Given P{K) satisfying Hypothesis \4-l[ CLnd given r > 0, 
there exists a right- continuous (for t > 0), time-homogeneous Markov pro- 
cess Xt with Xq = Xq such that 

supE'^'^le-''^ {K - Xry] = P{K) Vi^ > 

r 

and such that {e~'^^Xt)t>Q is a local martingale. 



11 



Remark. Although we wish to work in the standard framework with right- 
continuous processes, in some circumstances we have to allow for an im- 
mediate jump. We do this by making the process right-continuous, except 
perhaps at t = 0. At t = we allow a jump subject to the martingale 
condition E[Xo] = xq. 

To exclude a completely degenerate case, henceforth we assume P{xq) > 
0. If P{xq) = then necessarily, to preclude arbitrage, P{K) = {K — xq)"^ 
and Xt = xoe*"* is consistent with the prices P{K). In this case r = is an 
optimal stopping time for every K. 

Given P{K) satisfying Hypothesis 14.11 and such that P{xq) > 0, define 
ip by 

K — X 

ip{x) = sup (26) 

K:K>x P\J^) 

for X S (0,Xo]. By Lemma [T^ ip : {0,xq] [1, C)o] is a convex, decreasing, 
non- negative function with (/^(xo) = 1. Further, 

ip{xo) - ipixQ -e) <1- ^ ^o + e 



K* — Xq K* — Xq 

SO ip'{xQ) < —1/{K* — Xq) < 0. For some values of x, the supremum in 
()26p may be infinite since P may vanish on a non-empty interval {0,K]. We 
define 

X = mi{x > : ip{x) < oo}, 

and in the case where x > we see that (p{x) = oo for x < x. In fact x > 
if and only if K_> 0, and then it is easy to see that these two quantities are 
equal. 

We extend the definition of ip to (xq, oo) in any way which is consistent 
with the convexity, monotonicity and non-negativity properties and such 
that lim^^^oo ¥'(x) = 0. It is convenient to use ip{x) = (x/xq)"^ i^o-)^o ^ for 
then 95' (x) is continuous at xq, and (p is twice continuously differentiable and 
positive on (xo,cxd). 

Given cp, define s : (x, 00) 1— > (—00,00) via s'(x) = (p{x) — xip'{x), with 
s(xo) = 0. Then s is a concave, increasing function, which is continuous on 
(x, 00) . (It will turn out that s is the scale function, which explains the choice 
of label.) The function s has a well-defined inverse (7 : (s(x),s(oo)) — > (x, 00), 
and if s(x) > —00 then we extend the definition of g so that g{y) = x for 
y ^ s{x). Note that g : (— oo,s(oo)) (x, 00) is a convex, increasing 
function with g{0) = xq. Also define f{y) = ip{g{y)). Then / is decreasing 
and convex, with /(O) = v'(xo) = 1. 

Example. For geometric Brownian motion, provided /3 7^ 1 we have 
Kx) = x(^(xi-'3-x^^)(l + /?)/(l-/3), 
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y{i-f3) ' 

xo(l + /3)_ 



1 1/(1-/3) 



g{y) = xo 1 + 



and 



xo(l + /3) 



1 -/3/(l-/3) 



1 + 



If /9 = 1 the corresponding formulae are s{x) = 2xoln(x/xo), g{y) = 



Remark. Recall that a scale function is only determined up to a linear 
transformation. The choice s{xq) = is arbitrary but extremely convenient, 
as it allows us to start the process Z, defined below, at zero. The choice 
s'{x) = ip{x) — xip'{x) is simple, but a case could be made for the alternative 
normalisation s'{x) = {(p{x) — xif'{x))/{l — XQif'{xo)) for which s'(xo) = 1. 
Multiplying s by a constant has the effect of modifying the construction 
defined in the next section, but only by the introduction of a constant factor 
into the time-changes. It is easy to check that this leaves the final model Xt 
unchanged. 

Our goal is to construct a time-homogeneous process which is consistent 
with observed put prices, and such that e~'^^Xt is a (local) martingale. In 
the regular case we have seen how to construct a diffusion with these prop- 
erties. Now we have to allow for more general processes, perhaps processes 
which jump over intervals, or perhaps processes which have 'sticky' points. 
One very powerful construction method for time-homogeneous, martingale 
diffusions is via a time-change of Brownian motion, and it is this approach 
which we exploit. 

5 Constructing Time-homogeneous Processes as 
Time-changes of Brownian Motion 

In this section we extend the construction in Rogers and Williams [121 Sec- 
tion V.47] of martingale diffusions as time-changes of Brownian motion, see 
also Ito and McKean [8, Section 5.1]. The difference from the classical set- 
ting is that the processes defined below may have 'sticky' points and may 
jump over intervals. Since diffusions by definition are continuous, the result- 
ing processes are not diffusions, but one might think of them as 'extended 
diffusions', and they are 'as continuous as possible'. 

Let be a measure on R and let = {J-u)u>o be a filtration supporting 
a Brownian motion B started at with local time process L^. Define F to 
be the left-continuous increasing additive functional 



Xo& 



.y/i2xo) and f{y) 



= e 



2//(2xo) 




(27) 
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and let A be the right-continuous inverse to T, i.e. 

At = inf{u : r„ > t}. 

Note that F is a non-decreasing process, so that A is well defined, and At 
is an F-^-stopping time for each time t. Set Zq = and for t > set 
Zt = and J^t = J^At • ^o^e that Z is right-continuous except perhaps at 
t = 0. The process Zt is a time-changed Brownian motion adapted to the 
filtration F = {J't)t>o and subject to mild non-degeneracy conditions on v, 
(see Lemma l5 . 1 1 below) . the processes Zt and Z| — At are local martingales. 
Further, iiv{dy) = dy/j{y)'^ then F„ = £j{Br)~'^dr and At = j/^ -f{Zs)'^ds, 
so that Zt is a weak solution to dZt = ^{Zt)dWt, and Z is a diffusion in 
natural scale. The measure v is called the speed measure of Z, although, as 
pointed out by Rogers and Williams, u is large when Z moves slowly. 

The measure i' may have atoms, and it may have intervals on which it 
places no mass. If there is an atom at z then dT^u/du ^ whenever By^ — z, 
and then the time-changed process is 'sticky' there. Conversely, if u places 
no mass in (a, /5) then F is constant on any time-periods that B spends in 
this interval, and the inverse time-change A has a jump. In particular, Zt 
spends no time in this interval. If iy{{z}) = oo then F„ = oo for any u 
greater than the first hitting time by B of level z. Then A^o < so 
that if Z hits z, then z is absorbing for Z. The other possibility is that Z 
tends to this level without reaching it in finite time. 

Define Zi, S (0, oo] and z^ £ [— cxd,0) via 

Zu = mf{z > : J^((0, z]) = oo} and z^ = sup{z < : iy{[z, 0)) = oo}. 

The cases where 'z^ = or z^, = correspond to the degenerate case 
Xt = xqc^^ mentioned in the previous section, and we exclude them. The 
following lemma provides a guide to sufficient conditions for a time-change 
of Brownian motion to be a local martingale, and therefore provides insight 
into the constructions of local martingales via time-change that we develop 
in the next section. 

Lemma 5.1. Suppose that either 'z^^ < oo or v charges (a, oo) for each a, 
and further that either > —oo or v charges (— oo,a) for each a. Then 
Zt = BAt ^-^ ^ local martingale. 

Proof. See Appendix A. 2. □ 



6 Constructing the model 

We now show how to choose the measure v which gives the process we want. 
Define v via 
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and let y{dy) = oo for y < s(0) in the case when s(0) > — oo. Similarly, 
in the case when s(oo) < oo we set v{dy) = oo for y > s(oo). Where g' is 
absolutely continuous it follows that u has a density with respect to Lebesgue 
measure, but more generally (j28p can be interpreted in a distributional sense. 

Now, for this v we can use the construction of the previous section to 
give a process Zt. If we then set Xt = g{Zt), then, subject to the hypotheses 
of Lemma l5.ll Zt = s{Xt) is a local martingale, so that s is a scale function 
for X. The process X is our candidate process for which the associated put 
prices are given by P. 

Example. For geometric Brownian motion, 

Hdy) ^ ±9;^ ^ 1 P 

dy 2r g{y) 2rxl (1 + (5)^ 

In the case /3 = 1 this simplifies to 

v{dy) ^1^1 
dy S^qF 

where C = 2xq^. Then r„ = uC'^; At = tC^; Zt = Btc^ = CBt for a 
Brownian motion B and Xt = Xoe^*^"^^" = xqc^^^*. 

Recall that r„ = Jj^L^i'{dz) and let ^ be the first explosion time of F. 
Note that by construction F is continuous for t < and left continuous at 
t = ^. Since 1^ is infinite outside the interval [s(0),s(oo)] we also have the 
expression ^ = inf{n : Bu ^ [•s(O), s(oo)]} = ^^q) ^^m,oo)' '^^^ inverse scale 
function g is concave on (s(0), s(oo)), but may have a jump (from a finite to 
an infinite value) at s(oo). In that case we take it left-continuous at s(oo) 
so that we may have g := lim^^oo 9is{z)) is finite. 

For < n < ^, define M„ = e-''^^g{Bu) and iV„ = e-''^"/(5„). 

Lemma 6.1. M = (Mu)o<u<^ o-nd N = {Nu)o<u<( o'^e ¥^-local martin- 
gales. 

Sketch of proof: Suppose that (p is twice continuously differentiable with 
a positive second derivative. Then g is twice continuously differentiable. For 
u < ^, applying Ito's formula to = e~^^"g{Bu) gives 

e'^-dMu = g{Bu)dB^ + 

But, by definition, dVu/du = g" {Bu)/{2rg{Bu)), so M is a local martingale 
as required. 

A similar argument can be provided for the process N. For the general 
case, see the Appendix. □ 



xo(l + /3) 



— Z 



du. 
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Since M and N are non-negative local martingales on [0, ^) they converge 
almost surely to finite values, which we label and N^. In particular, if 
^ = H^^^y then = 0. However, if ^ = i/fj^^ then there are several 
cases. The fact that a non-negative local martingale converges means that 
we cannot have both < oo and g = lim^i-oo 9{s{z)) = oo. Instead, if 

< cxD then g < oo and = e~^'^^g. If = oo, and g < oo then = 0, 
whereas if = oo and g = oo then (Mu)„<^ typically has a non-trivial 
limit. Similar considerations apply to A'^. 

Recall that A is the right-continuous inverse to T and define the time- 
changed processes Mt = Ma^ and Nt = Na^ Note that these processes 
are adapted to F and that, at least for t < T^, we have TAt = t, Mt = 
e-''g{Zt) = e-^'Xt and Nt = e-''f{Zt) = e-'^V(^t)- 

If (s(0), s(oo)) = M then ^ = oo, = oo and M is defined for all t. 

If s(0) > — oo then we may have ^ = H^^^y In this case either = 

Tttb = oo, whence M is defined for all t as before, or Ft < oo. Then 

Mr^ = M^ = e~''^ig{B^) = 0, and we set Mj = for ah t > F^. It follows 
that Xt = for all t > F^ and is an absorbing state. 

Similarly, if s(oo) < oo then we may have ^ = H^^^y Then either 

F| = oo, whence M is defined for all t, or F^ < oo. In the latter case, if 
^ = H^^^^ < oo and Fg < oo then Mr^ = = e'^'^ig. We set Mt = 

for all t > F^, and it follows that for t > Fg, Xt := e^'^Mt = e'"(*-^e)g. Thus, 
for t > T^, X grows deterministically. An example of this situation is given 
in Example 18. 4l below. (In fact, the case where g < oo, which depends on the 
behaviour of the scale function s to the right of xq, can always be avoided 
by suitable choice of the extension to ip.) 

We want to show how M and X inherit properties from M. The key 
idea below is that, loosely speaking, a time-change of a martingale is again 
a martingale. Of course, to make this statement precise we need strong 
control on the time-change. (Without such control the resulting process can 
have arbitrary drift. Indeed, as Monroe [lOJ has shown, any semi-martingale 
can be constructed from Brownian motion via a time-change.) We have the 
following result the proof of which is given in the Appendix. 

Corollary 6.2. The process {e~'''^Xt)t>o is a local martingale. 

We can perform a similar analysis on and N and use similar ideas 
to ensure that is defined on R^. The proof that is a local martingale 
mirrors that of Corollary 16.21 

Corollary 6.3. The process {e~'^^(p{Xt))t>o is a local martingale. 
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7 Determining the put prices for the candidate 
process 

Recall the definitions of s, g and u via s'{x) = ip{x) — xip'{x), g = 
and i^{dy) = g" (y) / {2rg{y))dy . Suppose that Z is constructed from and 
a Brownian motion using the time-change F, and construct the candidate 
price process via Xt = g{Zt). By Corollary 16.21 the discounted price e~^^Xt 
is a (local) martingale. To complete the proof of Theorem 14.21 we need to 
show that for the candidate process Xt the function 

P{K) := supE[e-''^{K - 

T 

is such that P{K) = P{K) for all K > 0. 

Unlike the regular case the process X that we have constructed may 
have jumps. For this reason, for x < xq we modify the definition of the first 
hitting time so that = infju > : X^ < x}. 

Theorem 7.1. The perpetual put prices for X are given by P. 

Proof: Fix x S {x,xo). Suppose first that x is such that F is strictly 
increasing whenever the Brownian motion B takes the value s{x). Then 
Xh^ = X. More generally the same is true whenever — 6, s{x)]) > 

for every 5 > 0. By Corollary 16.31 we have that {e~'^^ip{Xt))t<Hx is a local 
martingale, and ip is bounded on [x, oo) so it follows that e-'^(*^^-)(^(XMHj 
is a bounded martingale and y?(xo) = E^o[e~^^^(/3(x)]. Hence, 

-vH^tr^ M /T^ ^V'(a^0) K-X 



P{K) > E^'o[e~'''^-(/^ - x)] = {K-x) 



ip{x) (fix) 



Otherwise, fix x^{x) = mi{w < x : z^((s(w), s(x)]) = 0} and x^(x) = 
supjtD > X : i'{[s{x), s{w)) = 0}. It must be the case that ip is linear on 
(x^(x), x'^(x)) and bounded on [x~(x),oo) and 



It follows that 



P{K) > max E^«[e"'^^-(i^-u;)] 

w£{x~ ,x+} 

K — w K — X 
= max — -— > — -— . 

w&[x-,x+} ip\w) (p[x) 

P{K) > sup = P{K). (29) 

x:x<xo '^{X) 

(Clearly, if x < x then {K — x)/ip{x) = 0, so the supremum cannot be 
attained for such an x.) 

To prove the reverse inequality, we first claim that the left derivative 
D~(p of the convex function (p satisfies 

D ip{xo) := hm = — . 30 

ej,0 e K* — Xq 
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To prove (fSOl) . first note that choosing K = K* in (1261) yields (p{x) > 
{K* -x)/{K* -xq), so 

D-^{x,) = lim - ^(^^ < 



x^xo Xq — x K* — Xq 

Conversely, note that for each 5 > there exists a non-empty interval {xq 
e, xo) on which 

K* -5-x 
^[.x) < — . 

K* — — Xq 

Consequently, for x G {xq — e, j;o) we have 

(/9(xo) - ip{x) ^ -1 



Xq — X K* — 6 — Xq 

Thus D~ip{xo) > —1/{K* — Xq) since > is arbritrary, so (f30]) follows. 
We next claim that for each fixed K < K* we have 

ip{x)>{K-x)^/P{K) (31) 

for all X. Clearly, this holds for x > K and for x < xq. Similarly, if 
Xq < X < K ^ then it follows from ()30p and the convexity of ip that 

K* -X K-x K-x 
' - K* -XQ - K -XQ - P{K) 

It follows from ()31|) and Corollary 16.31 that for any stopping rule r we have 

E^o[e-""(i^ - < P(i^)E^'»[e-''XXO] < P(K)(^(xo) = P{K). 

Hence P(i^) < P{K) for K < K* , and in view of ([29]), = P{K). 

For > iT* it follows from P{K*) = P{K*) = K* - xq, the convexity 
of P and Hypothesis 14.11 that P{K) = K — xq = P{K), which finishes the 
proof. □ 



8 Examples 

The following examples illustrate the construction of the previous sections. 
The list of examples is not intended to be exhaustive, but rather indicative 
of the types of behaviour that can arise. In each example we assume xq = 1. 



8.1 The smooth case 

We have studied the case of exponential Brownian motion throughout. It 
is very easy to generate other examples, for example by choosing a smooth 
decreasing convex function (with ip{xQ) = 1 and liuix-ioo V^i^) = 0)i ^^'^ 
defining other quantities from ip. 
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Example 8.1. Suppose (p{x) = (x + l)/(2x'^). Then, from ([3]) we obtain 

2 / \ 2x ~\~ 3 

a ix) = r X > 

^ ^ x + 3 

and, from ([8]) 

with P{K) = {K-l) for K > 5/3. 



8.2 Kinks in P 

If the first derivative of P is not continuous then we find that ip is hnear over 
an interval (a, /?) say. Then s' is constant on this interval and g is linear 
over the interval It follows that does not charge this interval, 

so that is constant whenever Bu G (s(a), s(/3)), and At has a jump. Then 
Zt jumps over the interval (s(a), and Xt spends no time in {a, (3). 

Example 8.2. Suppose P{K) satisfying Hvpothesis 14.11 is given by 

r K^/8 0<K < 27/32 

P{K) = < 4K^/27 27/32 <K<3/2 
[ {K-l) 3/2 <K. 

Then P is continuous, but P' has a jump at = 27/32. 
Using (j26|) we find that 



ip{x) 



2x-^ <x < 27/64 
X > 9/16. 



Over the region / = [27/64,9/16] ip is given by linear interpolation. The 
corresponding scale function is linear on /, and in the construction of Z, 
v assigns no mass to s{I). The process X is a generalised diffusion with 
diffusion coefficient given by (t{x) = \/2r for x < 27/64, o"(x) = y/r for 
X > 9/16 (strictly speaking there is some freedom in the choice of a for 
X > xq = 1, but the constant value y/r is a natural choice) and (t(x) = oo 
for X £ I. 



8.3 Linear parts to P 

In this case the derivative of (p{x) is discontinuous at a point 7 say. Then s' 
is also discontinuous at this point, and g' is discontinuous at 5(7). It follows 
that u has a point mass at 5(7), and that includes a multiple of the local 
time at 5(7). 
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Example 8.3. Suppose P{K) satisfying Hypothesis 14.11 is given by 
PiK) = 



8K^/27 <K <3/4 

{2K-l)/4 3/A<K<l 



1<K<2 
(K -1) 2<K. 

Then P is convex, but is hnear on the interval [3/4, 1]. 
Then 

x~-^/2 < X < 1/2 



X ^ X > 1/2 



where we have chosen to extend the definition of 93 to (1, 00) in the natural 
way. Then s{x) = 3 — 2 In 2 — 3/22; for x < 1/2 and s{x) = 21nx otherwise. 
It follows that g is everywhere convex but has a discontinuous first derivative 
at z = —2 In 2, and that the corresponding measure has a positive density 
with respect to Lebesgue measure and an atom of size r~^/12 at —2 In 2. 
In the terminology of stochastic processes the process Z is 'sticky' at this 
point — for a discussion of sticky Brownian motion see Amit [3j or, for the 
one-sided case Warren |14j . 

If P is piecewise linear (for example if P is obtained by linear inter- 
polation from a finite number of options) then ip is piecewise linear, s is 
piecewise linear, g is piecewise linear, and consists of a series of atoms. 
As a consequence the process Zt is a continuous-time Markov process on a 
countable state-space (at least whilst Zt < s{xq) = 0), in which transitions 
are to nearest neighbours only. Holding times in states are exponential and 
the jump probabilities are such that Zt is a martingale. 

In turn this means that Xt is a continuous-time Markov process on a 
countable set of points (at least whilst Xt < xq). 

Example 8.4. Suppose 

P{K) = 




This is consistent with a situation in which only two perpetual American 
put options trade, with strikes 1 and 3/2 and prices 1/3 and 2/3, in which 
case we may assume that we have extrapolated from the traded prices to a 
put pricing function P{K) which is consistent with the traded prices. 

Consider the process Xt with state space {0} U {1/2} U [2, co) and such 
that 

• at t = 0-|-, X jumps to 1/2 or 2 with probabilities 2/3 and 1/3, respectively; 

• if ever Xt^ > 2 then thereafter Xt = Xt^e^'^*-*''); 

• zero is an absorbing state for X; 
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• if ever X reaches 1/2 then it stays there for an exponential length of time, 
rate 3r, and jumps to 2 with probability 1/3 and zero with probability 2/3. 

Note that for the continuous time Markov process Xt, then conditional 
on Xt = 1/2, 



limlE[X,,.-X,] 



3r 



i(2-X,) + ^(-X,) 



- = rXt 
2 



Also, for this process 
P{K) 

= max{K -1,{2K -l)/3,K/3} 
so we recover the put price function given at the start of the example. 



max Ele-'^^iK-XrY 

T=0,ifi/2,-ffo 



8.4 Positive gradient of P at zero, i.e. -P'(O) > 0. 

In this case limxio ip{x) < oo. It follows that s(0) > — oo and the resulting 
diffusion Xt can hit zero in finite time. In this case we insist that is an 
absorbing endpoint for the diffusion. 

Example 8.5. Suppose that 

r K/2 K < 1 

P{K) = l {K + lf/8 1<K<3 
[ K-1 K>3. 

Then (p{ x) = 2{x + 1)'^ an d x^aix)"^ = r{x + l){2x + 1), so that dXt = 
rXt + y'r{Xt + l){2Xt + l) dBf 

The following example covers the case of mixed linear and smooth parts 
to P{K), and shows an example where reflection, local times and jumps all 
form part of the construction. 

Example 8.6. Suppose P{K) satisfies 

' K/4: K <1 

K'^/A 1<K<2 
(K-1) K>2. 

Note that P' has a jump at x = 1/2. We have (p{x) = 4 — 4x for a; < 1/2 
and (p{x) = 1/x for 1/2 < x < 1. We assume this formula applies on [l,oo) 
also. 

It follows that for x > 1/2, rj{x)'^ = {xa{x))'^ = 2rx^. Note that since 
<y9(l/2) < oo we have -ffi/2 is finite (almost surely). Hence 

dXt = rXtdt + V2^XtdBt t < (32) 
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is consistent with the observed put prices, but we need to describe what 
happens when X hits 1/2. Note that ([32]) has solution Xt = e^^* for 

t < Hi/2- 

The process is not absorbed at 1/2 as this would conflict with the idea 
that the discounted price process is a local martingale, or equivalently dXt = 
rXtdt + dAdt, for a local martingale Mj. However, adding a reflection alone 
would again destroy the martingale property of e~'^^Xt. What the general 
construction above does is that a local time reflection and a compensating 
downward jump are added. This jump takes the process to zero where it is 
absorbed. 

Alternatively, the process can be formalised as follows. Let It = — infu<j{(i? 
In 2/\/2r) AO}. Then by Skorokhod's Lemma, Bt + It is a reflected Brownian 
motion (reflected at the level -(ln2/^/2r)) and e^(-S*+-f*) > 1/2. 

Let N'^ be a Poisson Process rate A, independent of B, and let be 
the first event time. Then the compensated Poisson Process {N^ — Xt)t>o 

also a martingale. 

Take A = \/2r and define X via Xq = 1 and 



and the compensated Poisson Process stopped at the first jump (A''^^ 
X{t AT^))f>o are martingales. The time change (N^ rpx — A(/t AT'^))t>o, is 



dXt = rXtdt + V2^Xt f dBt + dit t:It< 

Note that at the first jump time of the time-changed Poisson process, X 
jumps from 1/2 to zero. 

By construction (e~'^^Xt)t>o is a martingale. 

8.5 P is zero on an interval: P{K) = for K < K 

Now we find that '^{x) = oo for x < x where x = K_. Depending on whether 
the right derivative P'{K_+) is zero or positive, we may have ^{x+) is infinite 
or finite. In the former case we have that Xt does not reach x in finite time. 
In the latter case Xt does hit x in finite time. 

The first example is typical of the case where ^{x+) = oo, or equivalently 
where there is smooth fit of P at K_. 

Example 8.7. Suppose Xq = 1 and that P{K) solves 

K < 1/2 

P{K) = { {2K-lf/8 1/2 < A" < 3/2 




K > 3/2. 

Then P' is continuous and for 1/2 < x < 1 we have <^(x) = (2x — 1)~^. We 
also have f{x) = oo for x < 1/2. It does not matter how (p is continued on 
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(1, oo) (provided it is decreasing and convex), but for definiteness we assume 
the formula ip{x) = (2x — applies there also. 

It follows that 7]{x)'^ = {xa{x))'^ = r{2x — l)(4a; — l)/4. Note that since 
93(1/2) = 00 we have (the first hitting time of 1/2) is infinite. Hence 



dX, = rX.dt + ^ !l^^Jl(l^rf5^ t < H^^^ 

is consistent with the observed put prices, and since the process never hits 
1/2, it is not necessary to describe the process beyond ^^1/2- 

Now consider the other case where P'(K) > 0. 

Example 8.8. Suppose Xq = 1 and that P{K) solves 



P{K) 



' i^" < 1/2 

{2K-l)/A 1/2<K<1 

1<K<2 
K-l K >2. 



Then P' has a jump at = 1/2. 

We have (p{x) = 00 for x < 1/2 and f{x) = 1/x for 1/2 < x < 1. We 
assume this formula applies on [l,oo) also. 

It follows that for x > 1/2, t/(x)^ = {xa{x))'^ = 2rx'^. As before we have 

dXt = rXtdt + V2^XtdBt t < Hy2 (33) 

is consistent with the observed put prices, but we need to describe what 
happens when X hits 1/2. 

The probability that the process ever goes below 1/2 is zero, else the 
put with strike K = 1/2 would have positive value. The process cannot 
be absorbed at 1/2 as this would conflict with the idea that the discounted 
price process is a martingale, or equivalently dXt = rXfdt + dMt, for a local 
martingale Mf. Instead, the time-change T includes a multiple of the local 
time at s(l/2) = -2 In 2. 

In fact g" /g is constant for y > s(l/2) and equal to zero for y < s(l/2). 
The process is a linear combination of and L"^'"^^, where is the 
amount of time spent by the Brownian motion above s(l/2) before time u. 
It is easy to check using Ito's formula that e~^^^g{Bu) is a martingale in 
this case. The process Zt = B^t is 'sticky' at s(l/2) (this time in the sense 
of a one-sided sticky Brownian motion, see Warren [TJ]) and this property 
is inherited by X = s{Z). 

There is a third case, where ^{x_+) < 00, but if'{x+) = 00. 
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Example 8.9. Suppose ip{x) = 2 — ^/2x — 1 for 1/2 < x < 5/2. Equiva- 
lently, 



Then for 1/2 < x < 5/2 we have 

r]{xf = {xa{x)f = 2r{2x - 1) [2y/2x -1 + 1 - x) 

It follows that although Xt can hit 1/2, the volatility at this level is zero, 
and the drift alone is sufficient to keep Xt > 1/2. 

8.6 Kink in P at K*: K* < oo and P\K*-) < 1. 

In this case (p'{x) is constant on an interval (x,xo). This case is just like 
Section 

9 Extensions 

9.1 No options exercised immediately 

In Hypothesis 14.11 in addition to (i) and (ii) which are enforcible by no- 
arbitrage considerations, we also assumed (iii) that there exists a finite strike 
K* such that for all strikes K > K* the put option is exercised immediately. 
Since K* < cxo is equivalent to ip'{xQ) < 0, it is apparent from the expression 
in Q, that provided a is finite for some x > xq, or equivalently u gives mass 
to some z > 0, then this property will hold. However, it is interesting to 
consider what happens when this fails. 

Suppose P{K) > K - xo for all K and lim^ P{K) - {K - xq) = 0. 
Then (^'(xq) = 0, but ip is strictly decreasing on {x,xo). The measure u 
places no mass on (0,oo), the process Zf spends no time on (0, oo) and Xt 
never takes values above xq. In particular Xt is reflected (downwards) at xq. 
The resulting model is consistent with observed option prices, but not with 
the assumption that the discounted price process is a martingale. However, 
by allowing non-zero dividend rates we can find a model for which the ex- 
dividend price process is a martingale, and for which the model prices are 
given by P{K). See Section [9^2] below. 

Now suppose limx P{K) — {K — xq) = 5 > 0. If P{xq) = xq then 
P{K) = K, and we have an extreme example which falls into this setting. 
For P as specified above we have that if{x) = 1 on (xq — 6, oo). The measure 
I' places no mass on {—S, oo) and s(x) = x — xq on this region. Except for 
time 0, the process Zt spends no time in {—6, oo) and Xt jumps instantly to 
Xo — 6, and thereafter spends no time above this point. 
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Note that if xq is not specified, then this case can be reduced to the 
previous case by assuming xq = K — Uuik P{K). This obviates the need for 
a jump at t = 0. 



9.2 Time homogeneous processes with non-constant interest 
rates and non-zero dividend processes 

In the main body of the paper we have assumed that the interest rate r is a 
given positive constant, that dividend rates are zero, and that o" is a function 
to be determined. However, the same ideas can be used to find other time- 
homogeneous models consistent with observed perpetual put prices, whereby 
the volatility function is given, and either a state-dependent dividend rate, 
or a state-dependent interest rate is inferred. 

Suppose X has dynamics dXt = {r{Xt) — q{Xt))Xt dt + Xtcr{Xt) dBf. 
Given put prices P{K) as before, define (p via ip{z) = miK:K>z{K—z)/ P{K). 
Then the relationship between ip and the characteristics of the price process 
X are such that ip solves L(p = 0, where L is given by 



Note that we now allow for any of cr, (7 or r to depend on x. To date we have 
assumed that r is constant and q is zero and solved for a, but alternatively 
we can assume cr{x) is a given function, and r is a positive constant, and 
solve for q, or assume that q and a are given, and solve for r. 

For example, if r and a are given constants then the (proportional) 
dividend rate process is given by 



If q is negative this should be thought of as a convenience yield. 

By allowing for dividend processes which are singular with respect to 
calender time, and which are instead related to the local time of X at level 
xq, it is possible to construct candidate price processes which spend no time 
above xq. For example, if L is the local time at 1 of X, and if 



then Xt reflects at 1, and if (^(x) = (E^[e"''^-])"^ for x < 1, then (p'{l-) = 0. 
This gives an example of a model consistent with the class of option prices 
described in Section 19.11 



Lu = -a{xfx^Uxx + x{r{x) - q{x))ux - r{x)u. 
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9.3 Recovering the model from perpetual calls 

The perpetual American call price function C : [0, oo) [0, xq] must be 
non-increasing and convex as a function of the strike K, and must satisfy 
the no-arbitrage bounds {xq — K)~^ < C{K) < xq- 

If there are no dividends (and if e^''"*Xt is a martingale) then the per- 
petual call prices are given by the trivial function C{K) = xq. 

So suppose instead that the (proportional) dividend rate q is positive. 
Let be the increasing positive solution to 

-a;^(7(x)^V" + — q)xip' — rip = 

normalised so that ■0(xo) = 1. Then, for z > x, W[e~'^^''] = '4){x)/'ip{z), and 
call prices in a model where dXt = (r — q)Xt dt + Xta{Xt) dBt are given by 

C'(ir) = sup E^o[e-"^(X^-i^)+] = sup ^"^r ■ 

r x:x>xo ij^i^) 

Example. Suppose X solves {dXf/Xt) = (r — q) dt + a dBf, with Xq = xq. 
Then V'(x) = (x/xq)'^ where 7 = /3+ and 




Note that since g > we have 7 > 1. Note also that (p{x) = x^-. 
The corresponding call prices are given by 

C{K) = x^ sup {{x - K)x-^} 

x:x>xo 

which for K < (7— l)a;o/7 gives C{K) = (xq — K), and for K > (7— l)a;o/7 

c{K) = x2;7-^(7 - ly-^K^-"'. 

The example discusses the forward problem, but the discussion of the 
inverse problem is similar to that in the put case. Given perpetual call prices 
C{K), for a; > a;o we can define ijj via ^(x) = ini K:K<x{x — K)/C{K), and 
then construct a triple (T{x),q{x),r{x) so that 

■^x'^a{x)'^ip" + {xr{x) — q{x))il)' — r{x)il) = 0. 

By combining information from put and call prices it is possible to determine 
a model which simultaneously matches both puts and calls. 
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A Proofs 
A.l Duality 



Proof of Lemma 12.51 It is clear that g is non-negative and non-decreasing 
since / is positive and non-increasing. The lower bound on g follows from 
choosing z = zq /\k \n ([TB]) . and the upper bound follows since / is non- 
increasing. To show that g is convex, first note that g{k) is minus the 
reciprocal of the slope of the tangent of the function / which passes through 
the point (A;,0). For two given points ki and k2 with ki < k2, let li{z) 
and hiz) be the corresponding tangent lines. Let k = Xki + {1 — X)k2 for 
some A € (0,1), and let l{z) be the line through the point (0, A;) and the 
intersection point of li and 12- If the intersection point is denoted {z,l{z)), 
then the convexity of / guarantees that 

k — z (1 — X)ki — (1 — X)z Xk2 — Xz , , /, N 

s -w = ^ — W — - + -hw- = - + 

which proves that g is convex. 

To prove the self-duality, let z < zq. By the definition of g we have that 
g{k) > {k — z)l f(z) for all k > z. Consequently, 

=sup^ </(z). 
k>z 9[k) 

For the reverse inequality, let z < zq and let I be a tangent line to / through 
the point {z, f{z)) (such a tangent is not necessarily unique if / has a kink 
at z). Assume that the point where / intersects the z-axis is given by (A;', 0). 
Then g{k') = {k' - z)/f{z), so 

7— 1 / \ ^ ^ ^ ^ r / \ 

F{z) = sup — — - > — — = /(z , 

k>z g{k) g{k') 

which finishes the proof of (ii). The proof of (iii) can be performed along 
the same lines. □ 



A. 2 Time changes of local martingales 

Proposition A.l. Suppose (7u)u>o is a martingale with respect to the filtra- 
tion G = {Gu)u>0! o^'^^d At is an increasing process such that At is a stopping 
time with respect to G for each t. Define ^t = lAn o.'i^d Qt = Gai ■ general 
{lt)t>o is not a martingale. However, if ^ is a hounded martingale then 7 
is a hounded martingale. 
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Proof. Given a Brownian motion B, for b > 0, let -ff^ be the first hitting 
time of level b. Then (-Bf,);,>o defined via = Bttb is not a martingale. 

b 

However, if 7 is bounded then E[7t|^s] = IE[7aJ^aJ = 7^^ = 7s by 
optional sampling. □ 

Suppose now that we are in the setting of Section [5] where Zt is con- 
structed from the Brownian motion B. In particular, is an increasing 
additive functional of B, and A is the right-continuous inverse to T. 

Proof of Lemma 15.11 Intuitively a time change of Brownian motion is a 
local martingale, but if the additive functional T is constant when B is in 
[a, 00) then the resulting process spends no time above a and reflects there. 
To maintain the local martingale property we need either that the time- 
changed process never gets to a, or that there are arbitrarily large values at 
which r is strictly increasing. 

If [z_i^,z,y] is a bounded interval, then A^o < ^Hf^ and {Zt)Q<t<oo = 
{BAt)Q<At<Aoo is a bounded martingale by Proposition lA. 11 

Now suppose {z_y.,'Zy) = M and suppose that for each a, u assigns mass 
to every set (a, 00) and (—00, a). 

We have A-p^ > t with equality provided T is strictly increasing at t. Let 
{a+} and {a~} be two sequences converging to +00 and —00, respectively, 
so that u assigns mass to any neighborhood of and o~, and set Hn = 
\ni{u : Bu ^ (a~,a^)}. Then F is strictly increasing at Hn- Set Tn = Fjj^. 
Then At^ = Hn- Note that F u increases to infinity almost surely, and hence 
Th^ I 00. Under our hypothesis, {Mf"-)t>o given by 

Mf" := Mmt„ = i?A„T„ = Ba,ah„ 

is a bounded martingale. Hence T„ is a localisation sequence for M. 

The mixed case can be treated similarly. □ 

Proof of Lemma 16.11 For y £ (s(0), s(oo)), set H{y) = ln{g{y) / xq) and 
write h{y) = H'{y) = g'{y)/g{y). Then 

'^(dy) = ^^dy = f {h'idy) + hiyfdy) 

and, as usual, v(dy) = 00 for ?/ ^ [s(0), s(oo)]. 
We have B{y) = h{v)dv = \n{g{y)/xo) and 

^n = ^ [ Ll{H'\dy)+H'{yfdy). 

Let ^ be the first explosion time of F. Then, by the Ito-Tanaka formula (eg. 
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Revuz and Yor \TT\ Theorem VI. 1.5]), for u < ^, 



HiBu) = r H'{B,)dB, + \ f H"{y)Lldy 
Jo ^ Jr 

h{Bs)dBs- \ f Llhiy fdy + rTu. 
) ^ Jr 



Thus g{Bu) = xoe^^^""^ = £{h{B) ■ B)ue'^^'^, where £ denotes the Doleans 
exponential, and e~^^"'g{Bu) is a local martingale. It follows that M is a 
local martingale. 

Now define J{y) = Jq j{v)dv = ln/(y) and 



f„ = l / Ll{j"{dy) + J'{yfdy) 
Jr 



By exactly the same argument as above we find that f{Bu) = e"^^^^^ = 
£{j{B) ■ i?)„e''^" and e~^^^f{Bu) is a local martingale. 

It remains to show that in fact r„ = F^. Define L{y) = {f{y)g{y))~^ ■ 
Then L'{y)/L{y) = -g'{y)/g{y) - f{y)/f{y) = -{H'{y) + J'{y)) and 



J'{y)-H'iy) 



v'{9{y))9'{y) g'jy) _ g' {yMy)v>' (giv)) - My))] 
v{g{y)) g{y) g{y)v{g{y)) 
g'{y)s'{g{y)) _ 



g{y)f{y) 

and then 

(J'(y) - H'{y))' = {H'{y) + J'{y)){L{y)) = H'{yf - J'{yf. 

Finally, since Lu is a bounded continuous function with compact support 
for each fixed u, we conclude that F^j = Ty^. 

Proof of Corollary 16.21 Recall that in our setting F defined via 
grows without bound and is continuous, at least until B hits s(0) or s(oo). 
Thus, if ^ denotes the first explosion time of F, then the inverse function A 
is defined for every t and At = for t > T^. Then, using the extension of 
the definition of M beyond F^ as necessary, we have 

^* - " ^* " I t > F^. 

Recall that f is extended to (xq, oo) in such a way that lim^^ioo V^i^) = 0- 
Therefore either s(oo) < oo and v assigns infinite mass to all points z > 
s{oo) = 'Zu, or s(cxd) = oo and there exists a sequence a,i j oo such that u 
assigns mass to any neighbourhood of a„. 
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Suppose, we are in the second case. If s(0) > — oo then H^,^^ = ^ < cc, 



else ^ = 00. On H^^ < -ff^g) = ^! -I^m is strictly increasing at u = H^^ and 
At.,, =H?. Set 



H. 



B 
an 



" " I > H^ior ^ ^ 

where the second line is redundant if s(0) = —00. Then At,^ = H^^ A ^ 
is such that M^^" := Mf^Tn = -^AtA^Ai/l* < dio-n) and T„ is a reducing 
sequence for M. 

Now suppose s(oo) < 00 and g = 00. Choose a„ ] s{oo) such that v 
assigns mass to any neighbourhood of a„. Then, on H^^^^ < -f^^g) have 
by the argument after Lemma 16.11 that Tjjb | 00 almost surely, and the 

tin 

argument proceeds as before with T„ given by (|34p a reducing sequence. 

Finally, suppose s{oo) < 00 and g < 00. Then M is bounded by g and 
M is a martingale. □ 
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